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ABSTRACT 

We study O(N) symmetric supersymmetric models in three dimensions at finite 
temperature. These models are known to have an interesting phase structures. In 
particular, in the limit N — > oo one finds spontaneous breaking of scale invariance 
with no explicit breaking. Supersymmetry is softly broken at finite temperature 
and the peculiar phase structure and properties seen at T = are studied here at 
finite temperature. 
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1 Introduction 

Bosons and fermions behave differently as they interact with an heat bath and 
supersymmetry is softly broken. Breakdown of supersymmetry at finite tempera- 
ture has attracted much attention since the early interest in supersymmetry and 
involved much controversy on its consequences, appearance and phase structure 
. Other related issues such as the restoration of broken internal symmetries 
at finite temperature supersymmetric theories were also debated @. More re- 
cently, there is a continued interest in the thermodynamics of supersymmetric 
Yang-Mills theory and temperature effects on the minimal supersymmetric 
model ||. 

Two supersymmetric models are studied in this paper at finite temperature, 
the O(N) symmetric supersymmetric $ 4 field theory, which is shown to have a 
peculiar phase structure in the large iV limit, and a supersymmetric non-linear 
cx-model. 

At zero temperature, it has been found in the past ||[F|,|8[ that the large iV limit 
of the scalar </> 6 theory in d = 3 dimensions and the supersymmetric $ 4 theory 
|§ in d = 3 posses spontaneous breaking of scale invariance unaccompanied by 
explicit breaking. Here, we find out which of the peculiar properties of the phase 
structure seen at T = in these models are maintained at finite temperature 
and how phase transitions occur as the temperature, rather than the coupling 
constant, is varied. 

2 Supersymmetric models in the large N limit 

A simple supersymmetric model, with one scalar superfield, in three euclidean 
dimensions at zero temperature is considered first. We then examine a super- 
symmetric non-linear a-model. 

2. 1 The supersymmetric scalar field in three dimensions 
We consider the O(N) invariant action 

S($) = Jd 3 xd 2 9 [iD$ • D$ + iW($ 2 /AT)] , (2.1) 

where $ is an iV-component vector: = f + Otp + ^66 F. Here, 6 is a 

two-component Majorana spinor, ip is an iV-component real scalar field, ip is 
an iV-component, two-component Majorana spinor and F is an iV-component 
auxiliary field. D = d/d6—ql 6 is the covariant derivative, the integration measure 
is d 2 9 = ^d02d9i and 9 a 9p = More details about the conventions used 

in the article can be found in the appendix. 
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In component form for a generic super-potential, one finds 

S = Jd 3 x [4* + i(y 2 -^'(^/iV)#-f/''(^/iV)(#)W)/iV 
+ i<p 2 U'\<p 2 /N)] • (2-2) 

The theory violates parity symmetry. Space reflection is equivalent to the change 
U^-U 

One of the subjects which will be of interest here is spontaneous breaking 
of scale invariance. When the binding potential U{<& 2 ) that binds bosonic and 
fermionic O(N) quanta is tuned to a particular strength at which O(N) singlet 
massless bound states are created, spontaneous breaking of scale invariance oc- 
curs and the resulting massless Goldstone particles appear as a supersymmetric 
multiplet of a dilaton and "dilatino". Though one may expect this as a result 
of non-renormalization of the coupling constant, in fact this happens also in the 
non-supersymmetric case. Explicit breaking of scale invariance will appear only 
at the next to leading order in 1/N as in the scalar case in the d = 3. At 
N — > oo flat directions in parameter space appear in the supersymmetric theory 
while maintaining the system's ground state energy at zero. 

2.2 Super -field formulation of large N steepest descent method. 

In the large N limit, a constraint on <E> 2 is introduced by adding to the initial 
action the term 

S L = J d 3 xd 2 6L(6) [$ 2 (6>) - NR(6)] , (2.3) 

and integrating over the superfields L(6) = M + 9£ + \69\ and R(6) = p + 
8a + \60s. N — 1 components of $ are first integrated out while keeping a 
test-component superfield <E»i = 0. One finds 

Z = J [d(j>] [dR] [dL] e~ SN (2.4) 

with the large iV action 

S N = J d 3 x d 2 9 [±D0D0 + NU(R) + L (</> 2 - NR)] 

+ i(AT- l)Str In [-DD + 2L] . (2.5) 

The two saddle point equations, obtained by varying the superfields <fi and R, 
are 

2L0-DD0 = O, (2.6a) 

L-U'(R) = 0. (2.66) 
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The third saddle point equation, obtained by varying L, involves the <p super- 
propagator given by the solution to the equation: 



(-DD + 2L(0)) A(fc, 0, 0') = o 2 (0' - 0), 



that is, 



A(fc,0,0) = ,o nT7^~T~\ + 13, I B<f 9 • l 2 - 7 J 



/c 2 + M 2 + A 
At coinciding arguments, one finds 

w) 1 + MSe 



k 2 + M 2 ' 



Mee 



k 2 + M 2 + X k 2 + M 2 



(2A 



The trace over the group indices yields a factor (N — 1). The third saddle point 
equation, obtained by varying L in Eq. ( |2.5[ ), is thus (for N ^ 1) 



R-cjP/N = — trA(fc,0,0) 



iV 



f2yr) 3 



d^A; 



1 + M00 



M00 



A; 2 + M 2 + A /c 2 + M 2 



(2.9) 



It is now convenient to introduce a notation for the boson mass: 

m ip = m = \/M 2 + A . (2.10) 
and write the saddle point equations in component form. Eq. ( |2.6a|) implies 



F - Mip = , 
Xip + MF = . 



(2.11a) 
(2.116) 



Eliminating F between the two equations one finds (pm 2 = and, thus, if the 
O(N) symmetry is broken the boson mass m<p vanishes. Eq. (|2.6fc|) yields 



M = U'(p), 
X = m 2 — M 2 = sU"(p). 



Eq. (|2.9|) in component form is 

p - <p 2 /N 1 
s - 2Fip/N 



d 3 p 



(2ir) 3 J p 2 + m 2 

2M r ,o / i 



(2tt) 3 



1 



p 2 + m 2 p 2 + 



(2.12a) 
(2.126) 



(2.13a) 
(2.136) 
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In these equations, a cut-off A is being introduced explicitly. We set 

' A d d k 



27r) d / k 2 + m 2 



n d (m). (2.14) 



Below, we need the first terms of the expansion of 0</(m) for m 2 — > 0. Introducing 
the cut-off dependent constant 



1 f A d 3 p 



rc (2tt) 3 J p 2 

we rewrite Eqs. Q2.13p as 



O 3 (0), (2.15) 



p _ ^/N = p c - -j-m , (2.16a) 

47T 

s-2F(p/N= — M(\M\ -m). (2.166) 



2tt 

Note that, here, the change U i— > —U corresponds to 

Fh-F, sh-s, Mh-M. 



Finally, we calculate the action density £ corresponding to the action Sn 
(Eq. ([2.5|)), £ = Sn /volume, for vanishing fermion fields: 



£/N = -\F 2 /N + \sU'{p) + MFlp/N + \\(p 2 /N - \Ms - \\p 

+ \ tr ln(-9 2 + M 2 + A) - \ tr ln(^ + M). (2.17) 

In d = 3, we have in Eq. ( |2.17| ): 

itrln(-a 2 + M 2 + A) - |trln(^ + M) = \p c \ - —— (m 3 - |M| 3 ) . (2.18) 

Note that the saddle point equations are recovered, in component form, by taking 
the derivatives of £ with respect to the various parameters. 

Using the saddle point Eqs. ( |2.6a| ), (|2.16| ) and ( |2.12a| ) to eliminate F,s,p, 
one eliminates also the explicit dependence on the super-potential U and the 
expression simplifies into 

£/N = \M 2 y 2 /N + -j-(m - \M\) 2 {m + 2\M\). (2.19) 

In this form we see that the action density £ is positive for all saddle points, 
and, as a function of m, has an absolute minimum at m = \M \, and thus A = 0, 
that is, for a supersymmetric ground state. 
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Eq. ( |2.6|) implies Mp = 0, and thus £ = 0. Therefore, if a supersymmetric 
solution exists, it will have the lowest possible ground state energy and any 
non-supersymmetric solution will have a higher energy. Since M, m, and tp are 
related by the saddle point equations, it remains to verify whether such a solution 
indeed exists. 

In the supersymmetric situation the saddle point equations reduce to s = 0, 
Mp = and 

p - Pc = p 2 /N - \M\/4n , M = U'{p). 

In the O(N) symmetric phase p = and \M\ = 4n(p c — p) = \U'(p)\. In the 
broken phase U'(p) = and p 2 /N = p — p c . We will show in the next section 
that these conditions can be realized by a quadratic function U(p) and, then, in 
both phases the ground state is supersymmetric and £ vanishes. 

2.3 The $ 4 super-potential in d = 3: phase structure 
We now consider the special case 

U(R) = fiR+ \uR 2 U'(R) = p + uR. (2.20) 

The dimensions of the 9 variables and the field $ are [9] = — \ , [<&] = 4 
[u] = . Power counting thus tells us that the model is renormalizable in three 
dimensions. Prior to a more refined analysis, one expects coupling constant and 
field renormalizations (with logarithmic divergences) and a mass renormalization 
with linear divergences. Using the solution for the two-point function W^ 2 \ in 
Eq. ( A1.6Q in the appendix, one infers that the coefficient A(p 2 ) has at most 



a logarithmic divergence, which corresponds to the field renormalization, while 
the coefficient C{p 2 ) can have a linear divergence which corresponds to the mass 
renormalization. 

For the quartic potential, Eqs. ( p.!2|) are now 



M = n + up, X = us. (2.21) 
We introduce the critical value of p: 

p c = -up c . (2.22) 

Taking into account Eqs. ( |2.21| ), one finds that Eqs. ( |2.1(]| ) can now be written 

as 

M = p — p c + up 2 /N — ^- a/M 2 + A , (2.23a) 

47T 

A = 2uMp 2 /N + —M (\M\ - \lM 2 + X) . (2.236) 



Eqs. ( |2.23D relate the fermion mass = \M\, the boson mass = \lM 2 + A 



and the classical field p. The phase structure of the model is then described 
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by the lowest energy solutions of these equations in the {fx — fj, c , u] plane. The 
invariance under the change U i— ► —U that was mentioned above and seen now 
in Eq. ( |2.20| ) and the equations that followed. This invariance is reflected into 
the phase structure of the model and one can restrict the discussion to u > 0. 
The phase structure is summarized in Fig. |I] in the first and second quadrant of 
the {[i - \i c , u/u c }. 

We find, indeed, that supersymmetry is left unbroken (A = 0) and the ground 
state energy £ = in each quadrant in the {/i — /i c , u] plane. This is consistent 
with Eqs. (|2.23| ) having a common solution with A = (thus = = \M\) 
and M<p = 0. They then reduce to 

M = fi — fj, c + Uj-jr — -^-|M| , Mip = 0. (2.24) 

The broken O(N) symmetry phase: The M = solution implies a sponta- 
neously broken O(N) symmetry, scalar and fermion O(N) quanta are massless 
and 

V 2 = -N{fi -n c )/u, (2.25) 

which implies that this solution exists only for \x < \x c . The solution exists in 
the fourth (and second) quadrant of the {fi — fi c ,u} plane. 

The O(N) symmetric phase: Here, (p 2 = and the equation 

M = /j,-fi c - (u/u c )\M\ (2.26) 

yields the common mass M for the fermions and bosons. In Eq. ( |2.26|) we have 
introduced the special value u c of the coupling u: 

u c = 4tt . (2.27) 

Eq. ( [2.26| ) splits into two equations, depending on the sign of M. The first 
solution 

M = M+ = (ji- n c )/(l + u/u c ) > , 

exists only for \i> \i c (first quadrant), as one would normally expect. 
The second solution 

M = M_ = (// - /i c )/(l - u/u c ) < , 

is very peculiar. There are two different situations depending on the position of 
u with respect to u c : 

(i) u > u c = 4-7T and then \i > \i c : the solution is degenerate with another 
O(N) symmetric solution M+. 
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Fig. 1 Summary of the phases of the model in the {fx — fx c ,u} plane. Here ra v = 
= \M±\ = ([Jb — li c )/(u/u c ± 1). The lines u = u c and \x — fi c = are lines of first 
and second order phase transitions. 

(ii) u < u c and then fx < fx c : the solution is degenerate with a solution of 
broken O(N) symmetry. 

The following different phases appear in the phase structure that is summarized 
in Fig. |l|: 

Region I : fx — fx c > , u/u c < 1: 
Here, there is only one 0(iV)-symmetric, supersymmetric ground state with 
= nritp = M + = (fx — fx c ) / (u/u c + 1) and ip 2 = 0. 

Region II : fx — fx c > , u/u c > 1: 
There are two degenerate O(N) -symmetric (<p 2 = 0), supersymmetric ground 
states with masses = = M + = (fx — fx c )/(u/u c + 1) and = = 
-M_ = (fi- fi c )/(u/u c -l). 

Region III : fx — fx c < , u/u c > 1: 
There is one supersymmetric ground state, it is an ordered state with broken 
O(N) symmetry (ip 2 ^ 0, = = 0). 

Region IV : fx — fx c < , u/u c < 1: 
There are two degenerate ground states: An O(N) -symmetric , supersymmetric 
ground states with masses = = M_ = (fx — fx c )/(u/u c — 1) and ip 2 = 0. 
The second ground state is a supersymmetric, broken O(N) symmetry state with 
= = and ip 2 ^ 0. 



The action density: In order to exhibit the phase structure in terms of the 
variation of the action density £, we plot in Fig. || the expression (|2.19|) , but 
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use only the fermion gap equation in Eq. ( |2.23a|) , in such way that £ remains a 
function of p and A, or equivalently p> and to = \J M 2 + A: 



1 1 

-£(m, if) = \M 2 (m, <p)^ + — [to - \M(m, p>)\} 2 {m v + 2 \M(m, <p)\) 



(2.28) 



0.15 




0.05 



Fig. 2 Region IV : The energy density W(m, p) 
a function of the boson mass (to) and A, where A 2 



^■£(m, p) as given in Eq. ( j2.28| ) as 



if /u c . Here \i — \i c = —1 , u/u c 



0.2. As seen here there are two distinct degenerate phases. One is an ordered phase 
(ip ^ 0) with a massless boson and fermion, the other is a symmetric phase (if = 0) 
with a massive (m = |M_|) boson and fermion. 



Fig. |] displays the action density £ restricted to (p = 0: 
11 2 

-rz£(m, ip = 0) = — — [to - \jj - ix c - [u/u c )m\] (to + 2 \ji - \i c - (u/u c )m\) . 

iV Z47T 

(2.29) 

Several peculiar phase transitions can be easily traced now in Eq. (|2.28|) . First, 
one notes the phase transitions that occur when /i — /i c changes sign. When 
< u < u c and \i > \i c , the system has a non-degenerate 0(iV)-symmetric 
ground state with bosons and fermions of mass M = M+. As — fi c changes 
sign (0 < u < u c fixed), two degenerate ground states appear, as seen in Fig. ^. 
Either M = and <p 2 = —N(fi — ii c )/u or the system stays in an O(N) -symmetric 
ground state with a mass \M- \ for the bosons and fermions. Similarly, when one 
goes from \i < fx c to \i > [i c at u > u c = 4n, the O(N) symmetry is restored but 
there are two degenerate ground states to choose from, M = M±. 
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Fig. 3 The action density W{m) = j^£{m,<p = 0) from Eq. ( 2,29| ) in region II of 



Fig. [T|. Here fi — fx c = 1 (sets the mass scale) and u/u c = varies between 1.6 and 
1.2. There are two massive, degenerate O(N) symmetric SUSY vacua with (p = and 
possibly \M-/M+\ >> 1 (as u — > u c ) with = m v = M+ = (fi — fj, c )/(u/u c + 1) 
and rriip = m v = — M_ = (li — fi c )/(u/u c — 1). 




Fig. 4 As in Fig.||, the action density W(m) = jj£(m, tp = 0) with fi — fi c = 1 and 
u/u c = 1.2 is changed to u/u c = 0.8 (from region II to region I of Fig. |l]). There are two 
degenerate O(N) symmetric SUSY vacua at u/u c = 1.2 with masses = m v = \M±\ 
where \M±\ = fi/(u/u c ± 1) while at u/u c = 0.8 there is a non-degenerate vacuum at 
= m v = M + . 



As seen in Figs. |3] and |], an unusual transition takes place when one varies the 
coupling constant u. For positive \i — [i c we find two degenerate ground states if 
u > u c . As u is lowered (at fixed fi — /i c ) the ground state with mass = = 
—M- disappears (\M-/M+\ — > oo) and only the O(N) -symmetric phase remains 
with M = M_|_. Namely, suppose we consider at {u > u c , fx — \i c > 0} a physical 
system in a state denoted by A and defined by {tp 2 = 0, M = M_}. Such a 
system will go into a state B defined by {ip 2 = 0, M = M_|_} when u decreases 
and passes the value u = u c = 4%. Now, if we consider the reversed process; 
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a physical system at {u < u c , (jl — [i c > 0} that is initially in the ground state 
B and u increases and passes u = u c . There is now no reason for the system 
to go through the reversed transition from B to A since the O(N) -symmetric 
states with M = M_ and M = M+ are degenerate and since supersymmetry 
is preserved the energy of state A will not decrease below zero. These peculiar 
phase transitions with (|M_/M + | — > oo) and with "infinite hysteresis" in the 
A — > B transitions are due to the fact that supersymmetry is left unbroken 
and the degeneracy of the two possible ground states in the leading order in 
1/N. If supersymmetry would have been broken by some small parameter, the 
lifted degeneracy would, most probably, have translated into a slow first order 
transitions between the otherwise degenerate ground states. 

The transition from the degenerate vacua at u/u c = 1.2 to a non-degenerate 
ground state at u/u c = 0.8 is shown in Fig. ^ (from region II to region I). This 
type of transitions between the different phases of Fig. [1] will be studied as a 
function of the temperature in the next section. 

Special situation. In general, when fj, = fj, c the mass M vanishes. However, 
there is a special case when 

u = u c = 4tv . 

Then the value of M is left undetermined. An accumulation point of coexisting 
degenerate ground states exist in the phase structure shown in Fig. [L| The case 
\i = n c represents a scale invariant theory where, however, the O(N) fermionic 
and bosonic quanta can have a non- vanishing mass m v = = \M\. Since u 
has not gone any perturbative renormalization there is no explicit breaking of 
scale invariance at this point. Thus, the only scale invariance breaking comes 
from the solution of the gap equation for M which leaves, however, its numerical 
value undetermined. We see a dimensional transmutation from the dimensionless 
coupling u that is fixed at a value of u = u c into an undetermined scale M. If 
M ^ the spontaneous breaking of scale invariance will require the appearance 
of a Goldstone boson at the point u = u c . The massless Goldstone boson is 
associated here with the spontaneous breaking of scale invariance |J. Moreover, 
since the ground state is super symmetric, we expect the appearance of a massless 
O(N) singlet Goldstone boson (a dilaton) and its massless fermionic partner 
(a "dilatino"). In order to verify these properties, we now calculate the (LL) 
propagator that will enable us to see these poles in the appropriate four-point 
functions. 

2.4 (LL) propagator, bound states and massless fermion and boson O(N) singlet 
bound states 

The (LL) propagator. We now calculate the (LL) propagator in the symmetric 
phase. Then, the propagators of the superfields <f> and L are decoupled. In 
the example of a quartic potential, the p field can be eliminated by gaussian 
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integration. The relevant part of the L-action then reads 
at r 

-— d 3 x d 2 9(L -fi) 2 + ±(N- l)Str In (-DD + 2L) . 
Zu j 

The calculation of the (LL) propagator involves the super- propagator (|2.7| ) . For 
the inverse propagator A^ 1 one finds 

a~ l \p) = — s 2 (o> -e)-2N j^A(k, 0, e')A( P - k, e, e>) 

with, here, (see Eq. ( |A1.6|) ) 

A(k, 9, 9>) = [1 + \M8\9> - 9)\ e^' . 

Then, 

[1 + M8 2 (9' -6)]e^ e ' 



A(k, 9, 8')A(p - k, 9, 9') 



(k 2 + M 2 )[(p + k) 2 + M 2 } ' 



Notice the cancellation of the factor e l9 ^ 9 which renders the integral more con- 
vergent that one could naively expect. The integral over k then yields the three- 
dimensional bubble diagram B(p): 

If d 3 k 1 

B{P) = Wfl (k 2 + M 2 )[(p + k) 2 + M 2 ] = ^ Ar ctan(p/2|M|) . 

At leading order for p small, we need only -B(O) = 1/8tt|M|. Then, 

A" 1 = --^ [1 + {M+ \M\(u c /u)}5 2 (9' - 6)] e^ e ' . 

Considering that in momentum space we have 

[-DD + 2fi}5 2 (9' -9) = 2 [-2 + 5 2 (9' - 6) ft] e~^ e \ (2.30) 

we conclude that for small M > the LL propagator corresponds to a super- 
particle of mass 2M(1 + u c /u). For M < the mass is 2|M(1 — u c /u)\. For 
\u — u c \ small, it is a bound state and at the special point u = u c the mass 
vanishes. 

More generally, we find 
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with 



m(p) = 2M + 



uB(p) 



We note that the only renormalization required at leading order is a mass renor- 
malization, a situation similar to the ip G scalar field theory . As a consequence 
field dimensions are not modified. 

Clearly, the propagation of the fields M(x) and £(x), as indicated in Eq. ( |2.31| ), 
when combined with the L(#)$ 2 interaction in Eq. ( |2.3|) , namely, 



J d 2 6L$ 2 = M{-\^ + Fip) - (piiP + iAv? 2 , 

describes the bound states in the 9999, ipip and tjj(p scattering amplitudes. For 
example ||, in the supersymmetric ground state case and with fx — fx c = 0, the 
-i/V scattering amplitude ^^ ^(p 2 ), in the limit p 2 — * satisfies 



2u 



1 + 



u 



M 



+ 



a 



4tt M 2n M 



4ni \M\ 



(2.32) 



for M < and w — > u c 

One notes here that the fermion massless bound state pole appears when a non- 
zero solution (M) to the gap equation exists (m v = = \M\) in the absence of 
any dimensional parameters (fi — \x c = 0) . This happens when the force between 
the massive ip and tp quanta is determined by u — > u c . The massless O(N) singlet 
fermion bound state excitation is associated with the spontaneous breaking of 
scale invariance. Similarly, the bosonic partner of this massless bound state 
excitation can be then seen, at the same value of the parameters, in ipup and ipip 
scattering amplitudes as Eq. ( |2.31|) shows. 

At \i = \i c in the generic situation M = 0, or for \p\ — > 00 we find 



B{p) 



1 



and, thus, 



N 



N 



A- L \p) = --5 2 (e'-e)-^-e* 



u 



4p 



As a consequence the canonical dimension of the field L is one, as in perturbation 
theory, and the interaction L$ 2 in Eq. (|2.3|) is renormalizable. 
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2.5 Dimensions 2 < d < 3 

In d space-time dimensions 2 < d < 3, the space of 7 matrices implies tr 1 = 2 
and one notices, after some calculations following the same pattern as discussed 
above, that Eqs. ( |2.13|) are replaced by 

P -tp 2 /N = n d (m), (2.33a) 
s - 2Fcp/N = 2M [fid(m) - n d (\M\)} . 

Taking into account the other saddle point equations and p c = O^(0), one obtains 
a generalization of Eq. ( |2.19| ): 

2S/N = M 2 cp 2 /N + (M 2 -m 2 )[n d (m) -Q d (0)] + J dX' Q d (\/ M 2 + X' 

(2.34) 

and, thus, 

S/N = \M 2 ^ 2 /N + K(d) [\{m 2 - M 2 )m d " 2 - (m d - \M\ d ) /d] , 

where 

r(i - d/2) 

K{d) = - (4tt)^ • 

Moreover, from Eq. (|2.34| ), 

has the sign of m — \M\ because Q' d (m) is negative. The function, therefore, has 
a minimum at m = \M\ for all values of d, and supersymmetry is maintained in 
the ground state for all 2 < d < 3. 

Two dimensions. It may be interesting to consider the same model in two di- 
mensions. The model now is super- renormalizable. Even at fi = it is not chiral 
invariant since a chiral transformation that changes "iJjiJj in —ifjifj, is equivalent to 
the change U 1— ► —U . 

The expression ( |2.35| ) is cut-off independent and has a limit for d = 2: 

S/N = \M 2 tp 2 jN + — \m 2 — M 2 — 2M 2 ln(m/M)l , (2.36) 

87T 

an expression which again has an absolute minimum at m = \M\. Then, a 
minimization with respect to M and <p yields M(p = 0. At the minimum £ 
vanishes. 

Taking into account A = and Eq. ( |2.33q| ) in the d = 2 limit, we obtain the 
gap equation 

M = n + uip 2 /N + — [ln(A/|M|] + O(l)) . 

Z7T 
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The solution M 2 = A = 0, therefore, is not acceptable. The O(N) symmetry 
is never broken. When the mass M is small in the cut-off scale, the l.h.s. is 
negligible and 

m = M oc Ae 2 "" /u . 
This solution exists only when —fi/u is positive and large. 

2.6 A supersymmetric non-linear a -model at large N 

We now consider the supersymmetric non-linear cr-model in d dimensions, 
2 < d < 3. The action 

= J- I d d xd 2 #D$ • D$ (2.37) 
2k J 

involves a scalar superfield <E>, an iV-component vector that satisfies 

$ • $ = N. 

This relation is implemented by introducing a superfield 

L(x, 9) = M(x) + 9£(x) + \ee\(x), 

where M(x), X(x) and £(x) are the Lagrange multiplier fields, and adding to the 
action 

s l = ^J d d x d 2 9 L(x, 9) [$ 2 (x, 9) - N] . (2.38) 
The partition function is given by (<S($, L) = S + Sl) 

[d$][dL] e~ s ^> L K 

In terms of component fields, the total action reads 

S=^J d d x{\^{-8 2 + M 2 )cp - + M)ip 

+ i X ( ( p 2 -N)-£$. (p )}. (2.39) 

As in the case of the $ 4 theory, we integrate out N — 1 superfield components 
leaving out = <&i: 

Z = J [d<j>][dL] e ~ SN ^ L \ 

where 

S N ((/), L) = - I d d x d 2 9 [±D0 • D0 + L ((f 2 - N)] + ^— ^Str ln(-DD + 2L). 
K J 2 

(2.40) 



16 



By varying the effective action with respect to 0, one finds the saddle point 
equation 

DD0-2L<£ = O, (2.41) 



which implies, for constant ip and if) = 0, 

F - Mip = and MF + \ip = . 



Varying the large iV action with respect to L and using the expression ( |2.7|) for 
the propagator, one finds (iV ^> 1) 



(2.42) 



k ~N 



d d k 
[2%)° 



1 + M99 



M66 



k 2 + M 2 + A k 2 + M 2 



(2.43) 



We now introduce the boson mass parameter m = V M 2 + A. In terms of its 
components Eq. ( |2.43| ) yields (for %j) = 0) 



1 - 



M 



N 



Kfld{m), 

M[fl d (\M\)-n d (m)] 



(2.44a) 
(2.446) 



Dimension d = 3. Introducing the critical (cut-off dependent) value k c : 

O 3 (0), 



1 



we can write Eqs. Q2.44 ) now as 



1 1 



K K c Nk 
,2 



m 
4n 



M 



1 



Nk 4tt 



—M{m- \M\). 



(2.45a) 
(2.456) 



The calculation of the ground state energy density of the non-linear a-model 
follows the similar steps as in the ($ 2 ) 2 model. One finds the action density 



S/N 



1 



2Nk 



2,„2 U m 2_ M 2 ) 



m (f 



1 1 



l 



12tt 



(m 3 - |M| 3 ) (2.46) 



and, taking into account the saddle point Eq. ( |2.44a| ) 



£/N = 2^ MV + " |M|)2(m + 2|M|) ' 



(2.47) 
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an expression identical to (|2.19|) , up to the normalization of (p. Again, if a 
supersymmetric solution exists it has the lowest energy. We thus look for super- 
symmetric solutions. 

In the O(N) symmetric phase (<£> = 0), Eq. (|2.456|) is satisfied while Eq. ( |2.45o| ) 
yields 

m=\M\= 4tc (— - - 

This phase exists for k > K c . 

In the broken phase m = M = 0, and Eq. ( |2.45fc|) is again satisfied. Eq. ( |2.45a| ) 
then yields 

ip 2 /N = 1 - k/k c , (2.48) 

which is the solution for k < K c . 

Since we have found solutions for all values of k, we conclude that the ground 
state is always supersymmetric. At k c , a phase transition occurs between a 
massless phase with broken O(N) symmetry (<^> 2 ^ 0) and a symmetric phase 
with massive fermions and bosons of equal mass. The phase structure of the 
supersymmetric non-linear a-model is analogous to the structure of the usual 
non-linear cr-model and less surprising than in the ($ 2 ) 2 field theory. 

Dimension d = 2. The phase structure of the supersymmetric non-linear 
cr-model in two dimensions is rather simple, and again analogous to the structure 
of the usual non-linear cr-model. Eq. (|2.45a|) immediately implies that M 2 +X = 
is not a solution and thus the O(N) symmetry remains unbroken (ip 2 = 0) for 
all values of the coupling constant k. Then, with a suitable normalization of the 
cut-off A, 

m = Ae- 2 ^ K . (2.49) 
Correspondingly, the action density becomes 

— = -— - — \(M 2 + A) ln(M 2 + A) - A In A 2 - M 2 In M 2 - A] (2.50) 

iV 2k 8tv 



and, taking into account Eq. ( |2.49|) 



£ 1 

— = — - [m 2 - M 2 + 2M 2 \n(M/m 

We recognize expression ( |2.36| ), and conclude in the same way that a supersym- 
metric solution has the lowest energy. Supersymmetry remains unbroken, and 
the common boson and fermion mass is 

M = m = Ae" 2 ^ K . (2.51) 



The condition McA implies that the non-trivial physics is concentrated near 
k = 0. 
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3 O(N) supersymmetric model at finite temperature 

We now study the phase structure of a supersymmetric model in d = 3 space- 
time dimensions at finite temperature. 

3. 1 The free energy at finite temperature 



We start again from the action ( |2.1| ), which is transformed into 



<S($, p,L) = J d 3 x d 2 9 [§D$ ■ D$ + NU($ 2 /N) + L{9){<$> 2 {9) - NR(6)] , 

$(6>, x), L(6, x) and R(9,x) being iV-component scalar superfields, parametrized 
as in the T = discussion. 

Integrating over N — 1 superfield components of $ and keeping $i = (the 
scalar component of the superfield <p is identified as f^i = <p) one finds an ex- 
pression for the partition function Z that differs from Eq. (|2.4j) only by the field 
boundary conditions 

Z= J[d(/)}[dR\[dL}e- SN ^' R > L \ (3.1) 

where Sn is the large N action fl2.5p: 

S N = J d 3 x d 2 9 [|D0D0 + NU(R) + L (<j> 2 - NR)] 

+ §(iV - l)Str In [-DD + 2L] . (3.2) 

Among the three saddle point equations at T = 0, only one is changed at finite 
temperature: 

2L(j> - DD</> = , (3.3a) 
L-U'(R)=0, (3.36) 

R-f/N = —txA(k,9,9). (3.3c) 
Eq. ( ggjj ) involves A(k,9,9) (Eq. Q) 

A ' fc - g ' = P + 4 + A + ^ ( P + M| + A - PTM?) ■ <X4) 

where we call Mt the expectation value of M{x) at finite temperature T. 
When written in components, Eq. ( |3.3a|) implies 

F - M T cp = , (3.5a) 
\<p + M T F = . (3.56) 
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and, thus, the Goldstone condition <p(X + Mj) = follows. 
Eq. ( p.36| ) implies 

M T = U'(p), X = sU"(p). (3.6) 

When calculating the trace in Eq. ( 13.3d, ), we have to take into account that 
bosons at finite temperature satisfy periodic and fermions anti-periodic boundary 
conditions. For bosons we introduce the function 

r ( t\= T V f A dd ~ lk 

A d^"^ 1/1 1 \ 

o + ZKraS ; • (3-7) 



(27r) d " 1 u(k) \2 eP u fW -1 



with u lfi (k) = a/ Mj, + A + k 2 . It is convenient here to introduce the boson ther- 
mal mass 



m T = ^M 2 +X. (3.8) 

The parameter tut characterizes the decay of correlation functions in space di- 
rections. Note that Mr, however, does not characterize the decay of fermion 
correlations, because fermions have no zero mode, the relevant parameter being 

y/M 2 +7T 2 T 2 . 

Combining Eq. (|3.7| ) with the 9 = part of the finite temperature Eq. ( |3.3dj 
we obtain 

2/*r f d 2 A; 1 (\ 1 



(27r) 2 ^(fc) V2 e^( fc )/ T -l 
= p c - ^ln(2sinh(m r /2T)), (3.9) 

where p c has been defined in Eq. (|2.15| ). 

In the same manner, for fermions we introduce the function 

02(M T ,T)= ———(-- . h . /rrt \ (3.10) 



27r) d - 1 u(k) \2 e ^W/ T +l 



with uj^(k) = ^Jk 2 + M|. (This is the fermion analogue of Eq. ( |3.7|) in which 
one recognizes again the sum of quantum and thermal contributions. Note that 
the function G2(Mj>,T) has a regular expansion in M| at M T = 0.) Combined 
with the 99 part of Eq. ( |3.3cf j , it leads to 

r d 2 k r i (\ i 



(2tt) 2 \2u) v (k) \2 e ^( fc )/ T -l 
1 1 



2aty(fc) V2 e^( fc )/ T +l 
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Integrating we obtain 

s - 2Fp/N = [ln(2 cosh(M T /2)) - ln(2 sinh(m T /2T))] . (3.11) 

The action density T at finite temperature, to leading order in 1/N, is given 
by(/3 = l/T) 

F = S N /V 2 p, (3.12) 

where is given by Eq. (|3.2|) at constant fields and V% is the two-dimensional 
volume. 

The supertrace term in Eq. (|3.2| ) decomposes into boson and fermion contri- 
butions. The boson contribution involves {u n = 2nirT): 

1 , , _o os r d^k 



Vr, 



d-1 



trln(-V 2 + m 2 T ) = J ^ ln(^ 2 + A; 2 + m|) 



ld-1 



where Vd-i is the (<i — 1) -dimensional volume. The result in Eq. (|3.13|) is given 
up to a mass independent infinite constant. 
Similarly, for fermions (now uj n = (2n + 1)ttT) 



Vd_itrl 



Thus, one finds 



trln(0 + M T ) = U JL* H<*i + k 2 + Ml) 
J ^ 1 r > nez 

I A^ln[2cosh(c^(fc)/2T)] . (3.14) 



— Str In [-DD + 2L] = — tr ln(-d 2 + M| + A) - — trln(^ + M T ) 
V2 V2 V2 



| ^ln[2sinh(/3^/2)]-2 | ln[2 cosh(^/2) 



ip c (m| - Ml) - (m 3 - |M T | 3 ) 
d 2 A; 



+ 2 / 7^9 I 111 ! 1 - e " /3 " ¥ '] - M 1 + e-^*]} . (3.15) 
J \^) 

The explicitly subtracted part reduces for Mt = M to the T = result. The 
other contributions to T are the same (up to the change M 1— > My) as in 
Eq. ( |2.17| ) and, therefore, one finds 

— ^ = -— + My — + A— + is(£7'(p) - M T ) - — (m£ - |M T | 3 ) 
N 2N N 2N 2 v J; 12tt V T 1 J 1 ; 

+ 1\( Pc - p )+tJ ^{l n [l-e-^]-ln[l + e-^]} , (3.16) 
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which reduces at T = to the action density of Eq. fl2.17j) , JF(T = 0) = £. 

Eq. ( p.!6| ) can be rewritten by using Eq. (|3.9|), which is also obtained by 
setting to zero the d/d\ derivative of Eq. ( |3.16| ), as well as Mt — U'(p) = from 
Eq. ( p.6|) (or equivalently, setting to zero the d/ds derivative of Eq. (|3.16|) ), one 
finds 



1 , 9 o? 2 1 

12 



T = \Mp^ + -(m T - \M T \) 2 {m T + 2|M T |) 



+ — ln(l-e-™ T / T )+T / —— {ln(l-e-^)-ln(l+e-^)}. (3.17) 

47T J (27TJ Z 

Eq. ( |3.17|) is the finite temperature version of Eq. ( |2.19|) . 

Inserting Eqs. Q and (gUD (with F = M T ip) into Eqs. Q with U(p) = 
pp + \up 2 one finds (p c = —up c ) 

2 

M T = p-p c + u^- ^-Tln(2sinh(im T /T)) , (3.18) 

which are the gap equations for T ^ 0. 

Solution to the saddle point equations: We first note that the r.h.s. of the 
gap Eq. (|3.18j ) diverges when — > 0. The finite temperature system in three 
dimensions has the property of a statistical system in two dimensions. In two 
dimensions a spontaneous breaking of a continuous symmetry is impossible, due 
to the IR behaviour of a system with potential massless Goldstone particles. 
Therefore, the O(N) symmetry is never broken, <p = 0, and thus, 

^ = ^ - ^ In (2 sinh (imr/T)) , (3.21)a 

M l~™J = I [In (2 sinh (Imr/T)) - In (2 cosh (±M T /T))] . (3.22)6 

Note that while the expression of T is complicated, its derivative with respect 
to remains simple: 

1 BF _ (ml - M|) 



N dm T lQnrriT tanh(mT / 2T) 

Therefore, the minimum still occurs at rriT = \Mj>\, but one verifies that tut = 
\Mt\ is not a solution to the saddle point equations. We have only found a lower 
bound on the free energy density J as a function of Mt'- 



T = NT 



J -^hitanh (\$sj M T + k^j . (3.24) 
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Its derivative with respect to Mt is 



-iV— — ^ lntanh(|M T |/2T). 



dM T 2n 

Therefore, the derivative has the sign of Mt- The lower bound has a limit for 
Mt — > 0, which is thus an absolute lower bound: 

^ = -^C(3)iVT 3 . (3.25) 

Similarly, the derivative of JF with respect to Mt at tut fixed is 
dT M T 



dM T 2?r 



[ln(2cosh(M T /2T)) - ln(2 sinh(m T /2T))] 



_m T - M T 
~ 2u ' 

We now have to find the solutions to the saddle point equations and compare 
their free energies. 

Note a first solution in the regime T — > 0, fi < [i c with \m,T\ <C T and 
\Mt -C T. Then, from (|3.1§| ) we find the boson thermal mass 

ttit ~ Te~ 2ir ^ c ~^/ uT . 

Eq. ( |3.19|) yields the fermion thermal mass 

M 7 



Since both rriT and Mt are very small, the free energy is very close to the 
lower bound ( |3.25[ ). The solution we have found here is a precursor of the zero- 
temperature phase transition, and corresponds to ip ^ in region III and IV of 
Fig. |. 

Other solutions exist in this regime but they converge, up to exponential cor- 
rections, to the finite masses of the T = spectrum and, thus, their free energy 
is much larger from the lower bound ( |3.24| ). 



It is expected that even for T larger and fi > fi c , the continuation of the small 
mass solution remains the solution with lowest energy. 



In the following we take m = tut as the free variable and, using Eq. (|3.18|) for 

Mt = My(m, (p, /i, u, T), we discuss 

W(m, ip, T) = —T{mT = m, M T {m, <p, /i, u, T), T). (3.26) 
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Since, as described above, at T ^ there is no breaking of the O(N) symmetry, 
we will discuss mainly (p = in Eq. ( |3.26| ). For T = 0, Eq. ( |3.26| ) results in 
jjFim, ip 2 , T = 0) = jj£(m,(p 2 ) which has been analyzed above. In particular 
we see the phase structure in Fig.|l] and the interesting degeneracy found in 
Regions II and IV and exhibited in Figs. |^-^. Rather than changing the value 
of [i — n c and the coupling u as has been done above at T = 0, we will be 
interested here to see the temperature dependence while these parameters will 
be held fixed. 

At finite temperature, the effective field theory describes the interactions of 
the fermions and bosons with the heat bath. This interaction acts like a source 
of soft breaking of supersymmetry. The short distance behavior is cured at finite 
temperature in a similar way it happens at T = 0. We now show that the general 
properties of the transitions as a function of the temperature T have a similar 
character as the transitions seen at T = when the coupling constant is varied. 

As seen in Fig. |, at T = there are two regions in the {fi — fi c , u] space where 
the vacua are degenerate. These are (region II : \i — \i c > and u > u c and 
region IV : [i — \i c < and u < u c ) . Fig. || shows the ground state free energy 
at T = in region IV. 

We discuss first region (II) fi — \i c > and u > u c : 
Here, as T increases the ground state with the smaller mass (ra^ = m = m + ) 
has a lower free energy than the heavier one (m v = m = m_ > m_|_) due to a 
sizable entropy contribution. (Fig. |B|). 

Peculiar transitions can occur in this system. If the system was initially at 
T = in the ground state with a boson mass ra^ = m_, it will eventually go, as 
T is increased, into the state with mass m v = m+, namely, into the lower mass 
ground state. This is due to the entropy negative contribution to the free energy 
forcing the system to favors the lowest mass state. On the other hand, a system 
that started at a high temperature in the state of low mass m v = m + will stay 
in this state as the system is cooled and will never roll back into the m = rri- 
high mass phase. Favoring the lowest mass phase as the temperature increases is 
a general effect that will occur in any physical system that is initially (at T = 0) 
mass degenerated. Here, supersymmetry imply that the m + and m_ vacua are 
at the same energy E = at T = 0. 

Region (IV) n < and u < u c : here one finds at T = two distinct degenerate 
phases. One is an ordered phase (if ^ 0) with a zero boson and fermion mass, 
the other is a symmetric phase (ip = 0) with a massive (m = m_) boson and 
fermion (as seen in Fig. [2]). 

As the temperature increases, the ordered phase with the broken O(N) sym- 
metry (m = 0, if 2 7^ 0) disappears into a ip 2 = symmetric phase and a very 
small mass ground state = m > is created (as seen in Fig. |6| and in Fig. [7]). 

Eventually, also the other O(N) symmetric vacuum (m = rri-,p 2 = 0) goes 
into the small mass 0(AQ-symmetric ground state. 

As in II, a system that was initially, at T = 0, in the 0(AQ-symmetric phase 
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Fig. 5 The ground state free energy -^!F{m = m v , ip, T) as given in Eq. ( 3.17| ) at p = 



as a function of the boson mass (m) at different temperatures. Here jA — fj, c = l(this 

sets the mass scale) , u/u c = 1.5 and T varies between T = 0.5. At T = two 

degenerate O(N) symmetric phases exist with a light m = m+ and heavier m = m— > 
m+ massive boson (and fermion). As the temperature increases, the light mass phase 
is stronger affected as its entropy increases faster and becomes the only ground state. 



(to. 



m. 



0) will decay into the smaller mass state when the temperature 



is raised. But upon cooling the system in the small mass phase will roll into the 
ordered phase (to^ = m = 0, ip 2 ^ 0) at T = 0. The system will never roll back 



into the symmetric phase (to v 



p = 0). 



4 The supersymmetric O(N) non-linear a-model at finite temperature 

The supersymmetric non-linear cr-model we consider here has already been dis- 
cussed above at zero temperature. The partition function of the O(N) super- 
symmetric non-linear cr-model in d dimensions is given by 



with 



Z = J [d$][dL] e - 5($ ' L) 

£($) = — [ d d x d 2 9 D$ ■ D$ + L($ 5 
2k J 



N) 



(4.1) 



(4.2) 



where k is a constant. The scalar superfield <&(x,6) in is an iV-component vector 
and the scalar superfield L(x, 9) is parameterized as in the T = case. After 
integrating out — 1 components of $, leaving out <j> = we obtain 



Z 



[d(f>}[dL] e 



-5at(<^,L) 
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Fig. 6 Same as Fig. ||| but the temperature has been increased from T = (in Fig. |2|) 
to T = 0.7 (here). j^!F(m = m v ,(p,T) as given in Eq. ( [3.171 ) as a function of the 
boson mass (m) and A, where A 2 = (p 2 /u c - Here fi — fx c = —1 and u/u c = 0.2. A 
non-degenerate O(N) symmetric ground state (99 = 0) appears with a very small boson 
mass (the non-zero mass is not seen here due to the limited resolution of the plot. See 
Fig. ). 

where 



S N ((f), L) = - I d d x d 2 6 [|D0 ■ D<P + L (<P 2 - N)] + — — ^Str ln(-DD + 2L). 
k J " 2 

(4.3) 

Differentiating the action Sn(<]>,L) with respect to the superfields (f> and L, one 
obtains Eq. ( |2.41| ) and the generalization of Eq. ( |2.43| ) : 



2L(f> - DD(p = 
N 6 2 



K 



N 



tiA(k,0,0). 



(4.4a) 
(4.46) 



The first equation is equivalent to F = Mt<-P and \<p + FMt = while the 
second equation can be compared with Eq. ( |3.3c| ), and thus leads to (for N 1 
and = 0) 



1 

k N 
2F(f 
~N~ 



d d k 1 (\ 1 

(2yr) d ZUk) \2 + e^( fc )/ T -l 



2M n 



d d k 
(2n 



(4.5a) 



1 



2u (p (k) 



1 



+ 



1 



» v (k)/T_i 
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Fig. 7 This figure displays the effect of increasing the temperature from T = in 
Fig. | to T that varies between T = - 0.7 (Fig. § has T = 0.7). The ground state 
free energy j^!F(m = m v , if, T) as given in Eq. ( |3.17[ )at ip = is plotted as a function 
of the boson mass (m) at different temperatures. Here /j, — fj, c = —1, u/u c = 0.2. At 
T = there are two degenerate phases: An 0(./V)-symmetric phase, shown here, with 
a massive (m = m_) boson and fermion and an ordered phase (tp ^ 0) with massless 
particles (both phases are shown in Fig. |2|). At finite temperatures the O(N) symmetry 
is restored (see Fig. ^) and a small mass ground state appears, the heavy mass state 
decays into the small mass ground state as seen here. 



1 



1 



3 ^(fc)/T 



(4.56) 



with Uyik) = y/rrij, + k 2 and u^(k) = a/ + k 2 . 



Dimension d = 3. As we have already discussed, the d = 3 finite temperature 
theory is analogous from the point of phase transitions to a two-dimensional 
theory. Therefore, the O(N) symmetry remains unbroken and tp = 0. We thus 
write the two gap equations only in this limit: 



1 



1 

K 



T 
2^ 

TM T 
2tt 



ln(2sinh(m T /2T)), (4.6a) 
[ln(2sinh(m T /2T)) - In (2 cosh(M T /2T))] . (4.66) 



The calculation of the ground state energy density to leading order for — > oo 
of the non-linear a-model at finite temperature follows steps similar to those of 
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the ($ 2 ) 2 field theory. The free energy at finite temperature is given by 



T = TS N /V 2 



(4.7) 



where Sn is given by Eq. (1~3) at constant fields and V2 is the two-dimensional 
volume. Here, the supertrace term in Eq. ( f4.^| ) has to be replaced by the finite 
temperature form as it appears in Eq. ( |3.16| ). The free energy then is given by 



+ T 



d 2 k 

(2tt) 2 



{Mi 



ln[l + e-^]}. 



(4.* 



The free energy in Eq. ( }4.8| ) can be compared with Eq. ( |2.46| ) which gives the 
ground state energy T at T = 0. 

It can be verified that, after using the zero temperature limit of Eq. ( |4.6a| ), 
the energy density becomes identical to the energy density of the ($ 2 ) 2 theory, 
up to a possible rescaling of the field. It is now clear that the same mechanism 
works here. Using Eq. ( |4.6a| ), one indeed finds the expression (|3.17| ) (for (p = 0): 



= i (mT _ i MT » 2 ( mT + 2 i Mt » + £ (m ^~ M ^ ln(1 - e " mT/T) 



+ T 



d 2 k 

(27T) 2 



{ln(l 



-/3uj v 



ln(l+e-^)}. 



(4.9) 



Solutions. The derivative of T with respect to \Mj\ at tut fixed is 
= \M^T [ln(2cosh(MT/2T)) _i n ( 2s inh(m T /2T))] 

It is convenient to introduce the notation 



X(k, T) = exp 



1 \K C K 



Then, the behaviour of T depends on the position of X with respect to 2: 

2tt / 1 1' 



X = 2 T 



In 2 V k 



K 



an equation which has a solution only for k > k c . 

For X < 2, the derivative vanishes at Mt = and is positive for all \Mj<\ > 0. 
Instead for X > 2, the derivative vanishes both at Mt = and 



\M T \ = 2Tln \UX + ^X 2 - 4) 



(4.10) 
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which are the two solutions of Eq. ( |4.6a|) . The minimum of T is located at the 
second value ( |4.10| ). 

We find, therefore, an interesting non-analytic behaviour: 

M T = for X < 2, 

M T = 2Tln [\{X + VX 2 - 4)] for X > 2 . 

Eq. ( [4.6a| ) then yields directly the boson thermal mass 



(X + ^X 2 + 4) 



rriT = 2T In 

For k < k c and T — > 0, we find the asymptotic behaviour 

rriT ~ TX(k, T), 

which is exponentially small, and Mt = 0. 

For k > k c and T — > 0, both mj and Mt converge toward the finite T = 
limit with exponential corrections. 

In the opposite high temperature limit T^> — l/« c |, we find that rriT is 
asymptotically proportional to T: 

m T ~2Tln((l + v / 5)/2), 

while M T = 0. 

It is not clear whether such a result can survive 1/N corrections. 

Dimension d = 2. In generic dimensions 2 < d < 3, a phase transition is 
not even possible at finite temperature and ip = (in two dimensions it is even 
impossible at zero temperature). The gap equations take the form 

- = n d (m T ) +T d - 2 f d (m 2 T /T 2 ), (4.11) 

K 

= M T [n d {\M T \) - T d - 2 g d (M T /T 2 ) - fi d (m) - T d ~ 2 f d (m 2 T /T 2 )] , (4.12) 
where f d and g d are defined by 

f x^ 2 dx 1 

f d (z) = N d - 1 7 == r~^= - 

Jo vr + ,2 exp|y ar + z| — 1 

= iVd " 1 /J (2/2 ~" )(d " 3)/2 ^ (413) 



and 



g d (s) = iV d _i 



00 x d " 2 dx 1 



o V x 2 + s exp 



oo 



tfd-i/ (^-«) (d - 3)/2 ^ I - (4-14) 
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For d = 2, it is then convenient to introduce the physical mass m solution of 

1/k = 2 (m), 

so that the equations become 

2 ( m ) = 2 (m T ) + f 2 (m 2 T /T 2 ), (4.15) 

and either M T = or 

n 2 (\M T \) - g 2 (M 2 /T 2 ) - 2 (m T ) - / 2 (m|/T 2 ) = . (4.16) 

The first equation which determines tut is identical to the equation one obtains 
for the usual non-linear cr-model. At low temperature my = m up to exponen- 
tially small corrections. At high temperature 

T 1 1 

~ — ln(mr/ra) ~ — In (T/m), 

HIT 7T 7T 

a consequence of the domination of the zero mode and the UV asymptotic free- 
dom of the non-linear a-model. 

Combining both gap equations we find 

n 2 (\M T \)-g 2 (M 2 /T 2 )=n 2 (m). 

The analysis of the equation for d = 2 shows that, for T large, the equation has 
no solution and thus = 0, but is has a solution for T small. The situation, 
therefore, is similar to what has been encountered in three dimensions. Again, 
an analysis of LL propagator and 1/N corrections is required to understand 
whether this result survives beyond the large N limit. 
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Appendices 

Al Conventions: Supersymmetry and Majorana spinors in d = 3 

We briefly recall the properties of Majorana spinors in three euclidean space 
dimensions and explain our notation. In three dimensions the spin group is 
SU{2). Then, a spinor transforms like 

%j )u = U%l}, UeSU{2) . 

The role of Dirac 7 matrices is played by the Pauli a matrices, 7^ = a^. More- 
over, 02 is antisymmetric while o"20> is symmetric. This implies 

0207*0-2 = - cr M U =a 2 Ua 2 ■ 

A Majorana spinor corresponds to a neutral fermion and has only two indepen- 
dent components ^1^2- 

The conjugated spinor is defined by ( T means transposed) 

i> = 7 V'02 <^> = ie a pil>0 , (A 1 - 1 ) 
(fa/3 = —£/3on £12 = 1) and thus ifj transforms like 

One also defines a spinor of Grassmann coordinates: 

6 = T 6a 2 , 

Since the only no n- vanishing product is 6*16*2, we have 

OaOp = \^ap0 ■ ■ 

The scalar product 6 ■ is 

9-9 = -216^2 Ojp = i<5 a/3 0i02 • 
If 6', 6' is another pair of coordinates, because o^o^ is symmetric, one finds 

Qa^O' = T 9a 2 a^' = , (Al.2) 
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and for the same reason 
Other useful identities are 

It is convenient to integrate over 9%, 62 with the measure 

d 2 9 = l -&6 2 &6 1 . 

Then, 

J d 2 99 a e p = ±d a ^ j d 2 ee-e = \. 

With this convention the identity kernel 5 2 (9' — 9) in 9 space is 

S 2 (9' -9) = (9' -9)-(9' -9). (A1.3) 

Super fields and covariant derivatives. A superfield &(8) has the expansion 

= + + \66F. (A1A) 

Again, although only two 9 variables are independent, we define the covariant 
derivatives D Q and D a (D = 02D): 

D 3 



" " 8S a 



D a = J- - (90) a . 



'a 



Then the anticommutation relation is 

{Da, D/?} = -2[0} a(3 . 

Also 

D Q D a = A-l - (90)a~ - 4r(0O) a + a 9 a d 2 
o9 a o9 n o9 n o9 a 



Since the cx M are traceless, using the identity (|A1.2| ) one verifies that DD can also 
be written as 

and, therefore, in component form 

D a D Q $ = 2F- 290ip + 99d 2 (p . 
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and 



(D Q D a ) 2 = 4d 2 



Super symmetry generators and WT identities. Supersymmetry is generated 
by the operators 

d9 a d9 a 
which anticommute with D a (and thus D a ). Then, 

{Q a ,Qf3} = 2l0Uf3- 

Supersymmetry implies WT identities for correlation functions. The n-point 
function W^ n \pk, 9k) of Fourier components satisfies 



Q a W [n) ; 
To solve this equation, we set 



k UUa 



W (n) (p, 9) = 0. 



W {n) (p, 9) = F {n \p, 6»)exp 



jk 



(A1.5) 



where is a symmetric function in the exchange {p^ 9i} <-> {pj, 9j}. It then 
satisfies 

£^ (n) (M) = °. 

k a 

that is, is translation invariant in 9 space. 

In the case of the two-point function, this leads to the general form 

W {2 \p, 9', 9) = A(p 2 ) [1 + C(p 2 )5 2 (9' - 9)} e 1 ^ 9 ' . (A1.6) 
= A(p 2 ) [1 + C{p 2 ){9' - 9){9' -9)+ %90 - \p 2 999'9'\ . 

Since the vertex functions satisfy the same WT identities, they take the 
same general form. 
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